DISCRETE TOMOGRAPHY OF ICOSAHEDRAL MODEL SETS 



CHRISTIAN HUCK 

Abstract. The discrete tomography of B-type and F-type icosahedral model sets is inves- 
tigated, with an emphasis on reconstruction and uniqueness problems. These are motivated 
by the request of materials science for the unique reconstruction of quasicrystalline struc- 
tures from a small number of images produced by quantitative high resolution transmission 
electron microscopy. 



1. Introduction 

Discrete tomography (the word "tomography" is derived from the Greek TO/ioa, meaning a 
slice) is concerned with the inverse problem of retrieving information about some finite object 
from (generally noisy) information about its slices. A typical example is the reconstruction of 
a finite point set in Euclidean 3-space from its line sums in a small number of directions. More 
precisely, a {discrete parallel) X-ray of a finite subset of Euclidean d-space R*^ in direction u 
gives the number of points of the set on each line in IR*^ parallel to u. This concept should 
not be confused with X-rays in diffraction theory, which provide rather different information 
on the underlying structure that is based on statistical pair correlations; compare [lO], [12] 
and [19]. In the classical setting, motivated by crystals, the positions to be determined form 
a subset of a common translate of the cubic lattice or, more generally, of an arbitrary 
lattice L in R^. In fact, many of the problems in discrete tomography have been studied 
on Z^, the classical planar setting of discrete tomography; see [21], [TT] and [16]. Beyond 
the case of perfect crystals, one has to take into account wider classes of sets, or at least 
significant deviations from the lattice structure. As an intermediate step between periodic and 
random (or amorphous) Delone sets, we consider systems of aperiodic order, more precisely, 
so-called model sets (or mathematical quasicrystals), which are commonly regarded as good 
mathematical models for quasicrystalline structures in nature |38j . 

Our interest in the discrete tomography of model sets is mainly motivated by the task of 
structure determination of quasicrystals, a new type of solids discovered 25 years ago; see [33] 
for the pioneering paper and [37l [25l [H] for background and applications. More precisely, 
we address the problem of uniquely reconstructing three-dimensional quasicrystals from their 
images under quantitative high resolution transmission electron microscopy (HRTEM) in a 
small number of directions. In fact, in [26] and [36] a technique is described, based on HRTEM, 
which can effectively measure the number of atoms lying on lines parallel to certain directions; 
it is called QUANTITEM (Quantitative ANalysis of The Information from Transmission 
Electron Microscopy). At present, the measurement of the number of atoms lying on a line 
can only be approximately achieved for some crystals; cf. [26j[36]. However, it is reasonable 
to expect that future developments in technology will improve this situation. 
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In this text, we consider both B-type and F-type icosahedral model sets A in 3-space which 
can be described in algebraic terms by using the icosian ring; cf. [8], [27] and [29]. Note 
that the terminology originates from the fact that the underlying Z-modules (to be explained 
in Section [3|) of B-type and F-type icosahedral model sets can be obtained as projections of 
body-centred and face-centred hypercubic lattices in 6-space, respectively. The F-type icosa- 
hedral phase is the most common among the icosahedral quasicrystals. Below, we nevertheless 
develop the theory for both the B-type (also called I-type) and the F-type phase. Well known 
examples of icosahedral quasicrystals include the aluminium alloys AlMn and AlCuFe; cf. [22] 
for further examples. 

In practice, only X-rays in yl-directions, i.e., directions parallel to non-zero elements of the 
difference set yl — yl of yl {i.e., the set of interpoint vectors of A) are reasonable. This is due 
to the fact that X-rays in non- yl-directions are meaningless since the resolution coming from 
such X-rays would not be good enough to allow a quantitative analysis - neighbouring lines 
are not sufficiently separated. In fact, in order to obtain applicable results, one even has to 
find yl-directions that guarantee HRTEM images of high resolution, i.e., yield dense lines in 
the corresponding quasicrystal yl. 

Any lattice L in R'^ can be sliced into lattices of dimension d — 1. More generally, model 
sets have a dimensional hierarchy, i.e., any model set in d dimensions can be sliced into model 
sets of dimension d — 1. In Proposition 13.161 it is shown that generic (to be explained in 
Section [3|) B-type and F-type icosahedral model sets can be sliced into (planar) cyclotomic 
model sets, whose discrete tomography we have studied earlier; cf. [Hill] and [23]. The latter 
observation will be crucial, since it enables us to use the results on the discrete tomography 
of cyclotomic model sets, slice by slice. 

Using the slicing of generic icosahedral model sets into cyclotomic model sets and the results 
from [4], it was shown in [24] that the algorithmic problem of reconstructing finite subsets 
of a large class of generic icosahedral model sets yl [i.e., those with polyhedral windows) 
given X-rays in two yl-directions can be solved in polynomial time in the real RAM-model 
of computation (Theorem 14. 3p . Since this reconstruction problem can possess rather different 
solutions, one is led to the investigation of the corresponding uniqueness problem, i.e., the 
(unique) determination of finite subsets of a fixed icosahedral model set yl by X-rays in a 
small number of suitably prescribed yl-directions. Here, a subset £ of the set of all finite 
subsets of a fixed icosahedral model set yl is said to be determined by the X-rays in a finite set 
U of directions if different sets F and F' in £ cannot have the same X-rays in the directions 
of U. Since, as demonstrated in Proposition 15.11 any fixed number of X-rays in yl-directions 
is insufficient to determine the entire class of finite subsets of a fixed icosahedral model set yl, 
it is necessary to impose some restriction in order to obtain positive uniqueness results. In 
Proposition 15.31 it is shown that the finite subsets F of cardinality less than or equal to some 
A; S N of a fixed icosahedral model set yl are determined by any set of A: + 1 X-rays in pairwise 
non-parallel yl-directions. Proposition 15.61 then shows that, for every R > and any fixed 
icosahedral model set yl, there are two non-parallel yl-directions such that the set of bounded 
subsets of yl with diameter less than R is determined by the X-rays in these directions. For 
our main result, we restrict the set of finite subsets of a fixed icosahedral model set yl by 
considering the class of convex subsets of yl. They are finite sets C C yl whose convex hulls 
contain no new points of yl, i.e., finite sets C C A with C = conv(C) Pi yl. By using the 
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slicing of generic icosahedral model sets into cyclotomic model sets again, it is shown that 
there are four pairwise non-parallel yl-directions such that the set of convex subsets of any 
icosahedral model set A are determined by their X-rays in these directions (Theorem I5.12|) . 
In fact, it turns out that one can choose four yl-directions which provide uniqueness and yield 
dense lines in icosahedral model sets, the latter making this result look promising in view 
of real applications (Example 15.141 and Remark IS.lSp . Finally, we demonstrate that, in an 
approximative sense, this result holds in a far more general (and relevant) situation, where 
one deals with a whole family of generic icosahedral model sets at the same time, rather than 
dealing with a single fixed icosahedral model set. 

2. Preliminaries and notation 

Natural numbers are always assumed to be positive, i.e., N = {1,2,3,...}. Throughout 
the text, we use the convention that the symbol C includes equality. We denote the norm 
in Euclidean d-space R,*^ by || • ||. The unit sphere in M!^ is denoted by S'^"^, i.e., S*^"^ = 
{x G ]R!^ \ \\x\\ = 1}. Moreover, the elements of S'^"^ are also called directions. Recall that a 
homothety h: — > IR°' is given hy x Xx + 1, where A G IR is positive and t G IR*^. We call 
a homothety expansive if A > 1. If x G IR, then [x\ denotes the greatest integer less than or 
equal to x. For r > and x G IR*^, Br{x) is the open ball of radius r about x. For a subset 
5 C R^, fc G N and > 0, we denote by card(5), J^{S), T<k{S), V<r{S), int(5), d{S), 
bd(S), conv(5), diam(5) and the cardinality, the set of finite subsets, the set of finite 
subsets of S having cardinality less than or equal to k, the set of subsets of S with diameter 
less than R, interior, closure, boundary, convex hull, diameter and characteristic function of 
S, respectively. The centroid (or centre of mass) of an element F G J-^SM^) is defined as 
{^f(zF f)/ card(F). A linear subspace T of R'^ is called an S-subspace if it is generated by 
elements of the difference set S — S := {s — s'\ s,s' G S} of S. A direction u G S*^"^ is called an 
S -direction if it is parallel to a non-zero element of S* — 5. As usual, denotes the group of 
units of a given ring R. Finally, for (o, b, cY G R^ \ {0}, we denote by the hyperplane 

in R"^ orthogonal to {a,b,cY. 

Definition 2.1. Let d G N and let F G J"(R'^). Furthermore, let u G S'^"^ be a direction 
and let be the set of lines in direction u in R'^. Then, the {discrete parallel) X-ray of F in 
direction u is the function X^F : £^ ^ No := N U {0}, defined by 

XuF{l) :=card(Fn^) = ^\p{x) . 

Moreover, the support {XuF)~^{¥l) of XuF, i.e., the set of lines in £^ which pass through at 
least one point of F, is denoted by s\vp\>{XuF). For z G R'^, we denote by 1^ the element of £^ 
which passes through z. Moreover, for 5 C R'', we denote by £f the subset of £^ consisting 
of all elements of the form where z G 5", i.e., lines in £^ which pass through at least one 
point of 5. 

Lemma 2.2. [IH Lemma 5.1 and Lemma 5.4] Let d G N and let u G S'^^^ he a direction. For 
all F,F' G r{W^), one has: 

(a) XuF = XuF' implies card(F) = card(F'). 

(b) If XuF = XuF' , the centroids of F and F' lie on the same line parallel to u. 
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Definition 2.3. Let d > 2, let [/ C S'^'^ be a finite set of pairwise non-parallel directions and 
let F G ^(R"'). We define the grid of F with respect to the X-rays in the directions of U as 

- n f u 

u&U \lesupp(XuF) j 
The following property follows immediately from the definition of grids. 

Lemma 2.4. Let d > 2. If U C S'^^^ is a finite set of pairwise non-parallel directions, then 
for all F,F' G J='{B.'^), one has 

{XuF = XuF' e U) ^ F,F' d Gfj = G^' . 

Definition 2.5. Let d > 2, let f C J^{W^), and let m G N. Further, let U C S"*"^ be a finite 
set of directions. We say that 8 is determined by the X-rays in the directions of U if, for all 
F,F' G (5, one has 

(X„F = XuF' yueU) =^ F = F' . 
Further, we say that £ is determined by m X-rays if there exists a set U of m pairwise 
non-parallel directions such that £ is determined by the X-rays in the directions of U. 

The following property is straight-forward. 

Lemma 2.6. Let d>2, let h: R"' — > R"' be a homothety, and let U C S'^"^ be a finite set of 
directions. Then, if F and F' are elements of T{'!R.'^) with the same X-rays in the directions 
of U , the images h{F) and h{F') also have the same X-rays in the directions of U . 

Gardner and Gritzmann introduced the so-called convex lattice sets, i.e., finite subsets C 
of some lattice L C R"^ with C = conv(C) H L; cf. [lH Section 2]. More generally, we define 
as follows. 

Definition 2.7. Let d G N and let S C R'^. A finite subset C of 5 is called a convex subset 
of S if it satisfies the equation C = conv(C) n S. Moreover, the set of all convex subsets of S 
is denoted by C{S). 

3. ICOSAHEDRAL MODEL SETS 

We shall always denote the golden ratio by r, i.e., r = {l + V5)/2. Mor cover, by .' we 
will denote the unique non-trivial Galois automorphism of the real quadratic number field 
(Q(r) = (Q(\/5) = (Q e Qr (determined hy y/b ^ -V^), whence r' = -1/r = 1 - r. Note 
that r is an algebraic integer (a root of X^ — X — 1 G Z[X]) of degree 2 over Q. Moreover, 
Z[r] = Z © is the ring of integers in Q(r) and, for its group of units, one further has 
Z[r]^ = {r" I s G 2} {i.e., r is a fundamental unit of Z[r]); cf [20]. 

3.1. Definition and properties of icosahedral model sets. Let H be the skew field of 
Hamiltonian quaternions, i.e., 

M = {a bi + cj -\- dk \ a, b,c,d £ R} , 

a four-dimensional vector space over R with a non-commutative multiplication determined 
by the following relations for the generating elements 1 (implicit in the above representation) 
and i,j, k: 

^ = j"^ = k^ = ijk = — 1 , 
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together with the requirement that IR is central in H. Note that IR, is precisely the center 
of H. The conjugate of a = a + bi + cj + dk £ H is defined hy a = a — bi — cj — dk, the 
reduced norm by nr(Q) = aa = a? + b'^ + (? + d'^ and the reduced trace by tr(a) = a + a = 2a. 
Moreover, we shall sometimes call Re(a) := a G IR the real part and Im(a) := {b,c,dY G IR^ 
the imaginary part of a. Let Hq be the set of quaternions with real part 0, i.e., 

Ho := {a G H I tr(a) = 0} = {bi + cj + dk\b,c,d £ B] ~ IR^ . 

The icosian ring I (c/. [8l[27l[29]) is the additive subgroup of H that is given by the integer 
linear combinations of the quaternions 

((±1, 0, 0, 0)*)'' , i ((±1, ±1, ±1, ±1)*)'' , i ((0, ±1, ±r', rYf , 

where we identify H with IR^ via the basis {l,i,j,k} and, as in [9l Chapter 8], the superscript 
A indicates that all even permutations of the coordinates are allowed. The members of I are 
called icosians. Note that I is a ring, because these generators (which have reduced norm 1) 
form a multiplicative group, the icosian group, of order 120. Note further that I is also a free 
Z[r]-module of rank 4. By [7], I is a maximal order of the quaternion algebra IH((Q(r)) over 
Q(r), defined similar to H as 

H((Q(r)) = {a + bi + cj + dk\a,b,c,d e Q(r)} . 

The set 

lo := Im(I n Ho) C IR^ 
of 'pure imaginary' icosians is generated as an additive group by the elements 

((±l,0,0)*)\i((±l,±r',±T)*)'' , 

where the superscript A is defined as above. Consider the standard body-centred icosahedral 
module AIb of quasicrystallography, defined as 

Mb := Z[r](2,0,0)*eZ[r](l,l,l)*©Z[r](r,0,l)* 

= Z[r](0,2,0)*©Z[T](-l,-T',r)*©Z[T](l,l,l)* 



/r/^^A^* At.'^G Z[r], with \. 
YPn,o) r'^i3 + rj + S = 0{mod2) j ' 



cf- |2l E] and references therein. One has Im(I) = ^A^b and, further, lo = ^Mp, where Aip 
is the standard face-centred icosahedral module of quasicrystallography, defined as 

/3,7,5 G Z[r], with 
/3 = r7 = (mod 2) 



Mf := (A 7, 5) 



(2) 



{(/3,7,5)* GTWbI P + J + 6= 0(mod2)} 
Z[r](2, 0, 0)* © Z[r](r + 1, r, 1)* © Z[r](0, 0, 2)* 
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= Z[r](0,2,0)*©Z[r](-l,-r',r)*©Z[r](2,0,0)* C Mb, 

where integers on top of the inclusion symbol denote the corresponding subgroup indices; 
cf- OE] again. Both Mb and Mp are free Z[r]-modules of rank 3, and are hence Z-modules 
of rank 6. Moreover, both A^b and Mf have icosahedral symmetry, i.e., they are invariant 
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under the action of the rotation group Y. This group is generated by the rotations which are 
given, with respect to the canonical basis, by the following matrices 




Note that Y is the rotation group of the regular icosahedron centred at the origin G 
with orientation such that each coordinate axis passes through the mid-point of an edge, thus 
coinciding with 2-fold axes of the icosahedron. Moreover, the matrix on the left (resp., on the 
right) is an order 2 (resp., order 5) rotation. 

Remark 3.1. There is another Z-module of rank 6, intermediate between M.p and Mb, 
which also has icosahedral symmetry. This is the standard primitive icosahedral module M.p, 
defined as 

A^p := {{(3,-f,6Y (EMbI P + j + 6= or T (mod2)} . 

In contrast to Alp and Mb, Mp fails to be a Z[r]-module. In fact, M.p is a Z[2r]-module 
only, and it is a Z-module of rank 6. 

By definition, model sets arise from so-called cut and project schemes; cf. [6l [27] for general 
background material and see [3] for a gentle introduction. In the case of Euclidean internal 
spaces, these are commutative diagrams of the following form, where vr and vrjnt denote the 
canonical projections; cf. [27] . 



^ — R'^ X R"^ ^ R'" 
(4) U U lattice U dense 

1-1 ~ 

Li — > L — > L* 

Here, L is a lattice in R"^ x R™. Further, we assume that the restriction 7r|^ is injective 
and that the image ■K\^t{L) is a dense subset of R"*. Letting L := vr(I/), the bijectivity of the 
(co-)restriction 7r|^ allows us to define a map .* : L — > R"* by a* := vrint((7r|^^) ^{a)). Then, 
one has L* = 7rint(-^^) and, further, L = {(/,/*) | / G L}. 

Definition 3.2. Given a subset C R"" with / mt{W) C W C cl(int(W^)) and 
cl(int(W)) compact, a so-called window, and any t G R"', we obtain a model set 

A{t,W) :=t + A{W) 

relative to the above cut and project scheme ((H) by setting 

A{W) := {a £ L\a* eW} . 

Moreover, R'^ (resp., R*") is called the physical (resp., internal) space. The map .* : L ^ R*", 
as defined above, is the so-called star map of A{t,W), W is referred to as the window of 
A{t,W) and L is the so-called underlying 'L-module of A{t,W). The model set A{t,W) is 
called generic if it satisfies bd(VF) r\ L* = 0. Moreover, it is called regular if the boundary 
bd(l^) has Lebesgue measure in R"*. 
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Remark 3.3. Every translate of a window W C is a window again. 

Definition 3.4. B-type icosahedral model sets A^^^{t,W) arise from the cut and project 
scheme dlj by setting d := m := 3, L := Im(I) and letting the star map .* : Im(I) be 
defined by applying the Galois conjugation .' to each coordinate of an element a G Im(I). We 
denote by the set of all B-type icosahedral model sets and define 2^ as the subset of all 
generic B-type icosahedral model sets. Additionally, for a window W C R^, we set 

lf{W) := {Af,,{t,s + W)\t,sen^} n 1^. 

F-type icosahedral model sets ^^o(^' ^) arise from the cut and project scheme (|4} by setting 
d := m := 3, L := Iq and letting the star map .* : Iq — > again be defined by applying 
the Galois conjugation .' to each coordinate of an element a G Iq. Moreover, the sets X^, Xj 
and Xj(Vr), where W CR-^ is a window, are defined similarly. Below, we say that Aicoit, W) 
is an icosahedral model set if Aico{t,W) = Af^^{t,W) or Aico{t,W) = Af^^{t,W). Finally, B- 
type (resp., F-type) icosahedral model sets are also referred to as icosahedral model sets with 
underlying TL-module Im(I) (resp., Iq). 

Remark 3.5. Both star maps as defined in Definition 13.41 are Q-linear monomorphism of 
Abelian groups and naturally extend to a monomorphism Q(t)^ — > R^, which we also denote 
by .*. Both in the B-type and the F-type case, we shall denote by .~* the inverse of the 
co-restriction of the corresponding star map .* : L — > L* to its image. The images of both 
maps 7: L ^ x R^, defined by a i-^ (a, a*), are indeed lattices in x R'^ ~ R^. In fact, 
these images have a natural interpretation as a weight lattice of type Z?g in the B-type case 
and a root lattice of type in the F-type case; cf. [U [9] for background. Finally, one can 
easily verify that the images Im(I)* and Iq* are indeed dense subsets of R^. 

We refer the reader to [27l[30] for details and related general settings, and to [6] for general 
background. Before we collect some properties of icosahedral model sets, recall the following 
notions. A subset A of R*^ is called uniformly discrete if there is a radius r > such that 
every ball Br{x) with x G R"' contains at most one point of A. Further, A is called relatively 
dense if there is a radius R > such that every ball Br{x) with x G R'^ contains at least one 
point of A. 

Remark 3.6. Let A be an icosahedral model set with window W . Then, /I is a Delone set in 
R^ (i.e., A is both uniformly discrete and relatively dense) and is oi finite local complexity (i.e., 
yl — yl is closed and discrete). Note that A is of finite local complexity if and only if for every 
r > there are, up to translation, only finitely many point sets (called patches of diameter r) 
of the form ylfl Br{x), where x G R'^; cf. [35l Proposition 2.3]. In fact, A is even a Meyer set^ 
i.e., yl is a Delone set and yl — yl is uniformly discrete; compare [27]. Further, A is an aperiodic 
model set, i.e., A has no translational symmetries. Moreover, if A is regular, A is pure point 
dijjractive, i.e., the Fourier transform of the autocorrelation density that arises by placing a 
delta peak (point mass) on each point of A looks purely point-like; cf. [35]. If A is generic, A 
is repetitive, i.e., given any patch of radius r, there is a radius R > such that any ball of 
radius R contains at least one translate of this patch; cf. [35]. If A is regular, the frequency 
of repetition of finite patches is well defined, i.e., for any patch of radius r, the number of 
occurrences of translates of this patch per unit volume in the ball Br{0) of radius r > about 
the origin approaches a non-negative limit as r ^ oo; cf. [3l]. Moreover, if A is both generic 
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Figure 1. A few slices of a patch of the icosahedral model set Af^^ (left) and 
their .*-images inside the icosahedral window in the internal space (right), both 
seen from the positive x-axis. 

and regular, and, if a suitable translate of the window W has full icosahedral symmetry (i.e., 
if a suitable translate of the window W is invariant under the action of the group Yj!^ of order 
120, where := Y* U {—Y*) and Y* is the group of rotations of order 60 generated by the 
two matrices that arise from the two matrices in ([3]) by applying the conjugation .' to each 
entry), then A has full icosahedral symmetry Yh := Y U {—Y) in the sense of symmetries of 
Ll-classes, meaning that a discrete structure has a certain symmetry if the original and the 
transformed structure are locally indistinguishable (LI) (i.e., up to translation, every finite 
patch in A also appears in any of the other elements of its Ll-class and vice versa); see [3] 
for details. Typical examples are balls and suitably oriented versions of the icosahedron, the 
dodecahedron, the rhombic triacontahedron (the latter also known as Kepler's body) and its 
dual, the icosidodecahedron. 

Example 3.7. For a generic regular icosahedral model set with full icosahedral symmetry Yh, 
consider Af^^ := 71^^(0,5 + W), where s := 10~^(1,1,1)* and W is the regular icosahedron 
with vertex set Yj^{t',0, 1)*; see Figure [T]for an illustration. 

3.2. Cyclotomic model sets as planar sections of icosahedral model sets. In this 
section, we shall demonstrate that both B-type and F-type icosahedral model sets A can be 
nicely sliced into cyclotomic model sets with underlying Z-module ^[^5], where the slices are 
intersections of A with translates of the hyperplane H^'^'^'^^ in orthogonal to (r, 0, 1)*. 
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Prom now on, we always let Cs := e^'^*/^, as a specific choice of a primitive 5th root of unity 
in C Occasionally, we identify C with R^. 

Remark 3.8. It is well known that the 5th cyclotomic field ^(Cs) is an algebraic number field 
of degree 4 over (Q. Moreover, the field extension Q(C5)/Q is a Galois extension with Abelian 
Galois group G((Q(C5)/(Q) ~ (Z/5Z)^, where a (mod 5) corresponds to the automorphism 
given by Cs Cf; c/. [HH Theorem 2.5]. Note that, restricted to the quadratic field Q(r), 
both the Galois automorphism of Q(C5)/Q that is given by Cs ^ Cf and its complex conjugate 
automorphism [i.e., the automorphism given by ^ Cf) induce the unique non-trivial Galois 
automorphism .' of Q(t)/(Q (determined by r i— > 1 — r). Further, ^[^5] is the ring of integers in 
^(Cs); cf. [311 Theorem 2.6]. The ring ^[Cs] also is a Z[T]-module of rank two. More precisely, 
one has the equality ^[Cs] = Z[r] © ^[rjCs; cf. [H Lemma 1(a)]. Since is also a primitive 
5th root of unity in C, one further has the equality ^[Cs] = ^[(5] = Z[r] © Z[r]C|. 

Definition 3.9. Cyclotomic model sets with underlying Z-module ^[^5] Acyc{t,W) arise from 
the cut and project scheme ([4} by setting d := m := 2, L := ^[Cs] and letting the star map 
: L ^ be either given by the non-trivial Galois automorphism of (Q(C5)/Q, defined by 
Cs ^ cf) or its complex conjugate automorphism. 

Remark 3.10. The star map as defined in Definition 13.91 is a monomorphism of Abelian 
groups. Further, the image of the map 7 : L ^ x R^, defined by a 1— > (a, a*^), is indeed 
a lattice in R^ x R^. Finally, one can verify that the image L*'^ is indeed a dense subset of 
R^. For the general setting, we refer the reader to [UIMIES]. By [241 Lemma 1.84(a)] (see 
also [23l Lemma 25(a)]), for all cyclotomic model sets A with underlying Z-module ^[(^5], the 
set of yl-directions is precisely the set of Z[C5]-directions. 

Example 3.11. For illustrations of cyclotomic model sets with underlying Z-module ^[(^5], 
see Figure [2] on the left and Figure [31 cf. Proposition 13.161 and Example 13 . 1 71 b elow . 

Lemma 3.12. For L £ {Im(I),Io}, the following equations hold: 

(a) Ln//("'0'i) = Z[r](0,l,0)*©Z[r]i(-l,-r',r)*. 

(b) (Lni?(^'0'i))* = n if(^'Ai). 

Proof. Part (a) follows from Equations ^ and ^ together with the relations Im(I) = ^M.b 
and lo = ^Mp. Part (b) follows from the identity ((r, 0, 1)*)* = (r', 0, 1)*. □ 

Definition 3.13. We denote by $ the R-linear isomorphism $ : ij('^'O'i) — > (C, determined 
by (0,1,0)* ^ 1 and ^( — 1,— T',r)* 1— > ^5. Further, will denote the R-linear isomorphism 
$*: ^ C, determined by (0,1,0)* ^ 1 and i(-l,-r,r')* ^ Cf- 

Lemma 3.14. The maps $ and are isometrics of Euclidean vector spaces, where H^'^''^'^\ 
and C are regarded as two-dimensional Euclidean vector spaces in the canonical way. 
Moreover, identifying C with the xy-plane in R'^, $ and extend uniquely to direct rigid 
motions o/R"^, i.e., elements of the group S0(3, R). 

Proof. The first assertion follows from the following identities: 

||r(0,l,0)* -hs^(-l,-T',r)*|| =\r + sC5\ = Vr"^ + - rsr' , 

||r(0, 1,0)* -Fsi(-1,-T,r')*|| =\r + sCi\ = V^r2 + 3"^- rsr . 
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Figure 2. The central slice of the patch of yl^^ from Figure [T] (left) and its 
.Mmage inside the (marked) decagon {s + W)riH^'^ '^'^^ (right), both seen from 
perpendicular viewpoints. 

The additional statement is immediate. □ 

Lemma 3.15. Let L G {Im(I),Io}. Via restriction, the maps $ and induce isomorphisms 
of rank two 'L[T]-modules: 

Proof. This follows immediately from the definition of <^ and <l>* together with Lemma 13.121 
and Remark 13.81 □ 



Proposition 3.16. Let A he a generic icosahedral model set with underlying li-module L, say 
A = A^co{t,W) . Then, for every \ £ A, one has the identity 

^[{A n (A + i?("'°'i))) -A) = {z G Z[C5] I G VFa} , 

where is the Galois automorphism o/(Q(C5)/Q, defined by ^5 ^ Cf (^f^d, 

Wx := ^*{{W n ((A - ty + i7("''0'^))) - (A - ty) . 

Thus, the sets of the form 

(5) -^{{A n (A + i/(^'°'^)))-A), 

where X G A, are cyclotomic model sets with underlying Ti-module ^[^5]. 

Proof First, consider ^>(^), where £ {An (A + F(^'°'^))) - A. It follows that fi G Lni?(^'°'^) 
and (// + (A - t)y = fj* + {X - ty £ W. Lemma [3T5l implies that ^>(^) G ^[Cs], say 
^(/i) = Q + (3(5 for suitable a,/3 G Z[r]. One has 

^^y^ = a' + p'cl = $*(//'^) G Wx . 
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Figure 3. Another two slices of the patch of Af^^ from Figured! 

Conversely, suppose that z £ ^[Cs] satisfies z*^ £ Wx. Then, there are suitable a,/3 G Z[r] 
such that z = a + j3C,^ and, consequently, z*^ = a' + f3'Ci ^ ^A- By definition of W\, one 
has z*^ = where ^ £ H^'^ satisfies /i + (A - t)* e W . Clearly, there exist r, s G IR 

such that ^ = r(0, 1,0)* + 1, —t,t'Y, whence $*(/^) = r + sCf - The linear independence 
of 1 and d over R now implies that r = a and s = /3, so that fi £ L*. Moreover, one can 
verify that one has £ (yln (A + - A and = a + (3C^ = z. This proves the 

claimed identity. The assertion is now immediate. □ 

Example 3.17. For an illustration of the content of Proposition [3TT61 in case of the icosahedral 
model set A^^^ from Example 13.71 see Figures [2] and [3l 

3.3. The translation module of icosahedral model sets. In order to shed some light on 
the set of yl-directions of an icosahedral model set A with underlying Z-module L, we first 
have to establish a relation between icosahedral model sets and their underlying Z-modules. 
We denote by m,- the Z[r]-module endomorphism of Q(t)^, given by multiplication by r, i.e., 
a I— > ra. Furthermore, we denote by m-j-* the Z[r]-module endomorphism of (Q(t)'^)*, given 
by a* 1-^ (ra)*. 

Lemma 3.18. The map mr* is contractive with contraction constant 1/r £ (0,1), i.e., the 
equality ||mT-*(a*)|| = (l/r) ||q;*|| holds for all a £ Q(t)^. 

Proof. For a £ Q(r)3, observe that ||m^*(a*)|| = ||(ra)*|| = ||r'a*|| = (l/r) ||a*||. □ 

Lemma 3.19. Let A be an icosahedral model set with underlying Z-module L, say A = 
^ico(i;W^)- Then, for any F £ J-{L), there is an expansive homothety h : R'^ — > R'^ such 
that h{F) C A. 

Proof. From int(W^) / and the denseness of L* in R'^, one gets the existence of a suitable 
ao £ L with ao* £ int(iy). Consider the open neighbourhood V := int(Vl^) — ao* of in 
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R^. Since the map m^* is contractive by Lemma fS. 181 (in the sense which was made precise in 
that lemma), the existence of a suitable A: G N is implied such that [nir*)^ {F*) C V . Hence, 
one has {{r^a + ao)* | a G F} C int(VF) C W and, further, h{F) C A, where h: ^ is 
the expansive homothety given by x i— > t^x + (ag +t). □ 

As an easy application of Lemma 13.191 one obtains the following result on the set of A- 
directions for icosahedral model sets A. 

Proposition 3.20. Let A he an icosahedral model set with underlying 'E-module L. Then, the 
set of A-directions is precisely the set of L-directions. 

Proof. Since one has A — A C L, every yl-direction is an L-direction. For the converse, let 
« G §^ be an L-direction, say parallel to a G L \ {0}. By Lemma [3.191 there is a homothety 
h: IR3 ^ j^s g^jpj^ ^i^g^^ h{{0,a}) C A. It follows that h{a) - /i(0) G {A - A) \ {0}. Since 
h{a) — h{0) is parallel to a, the assertion follows. □ 



4. Complexity 

In the practice of quantitative HRTEM, the determination of the rotational orientation of a 
quasicrystalline probe in an electron microscope can rather easily be achieved in the diffraction 
mode. This is due to the icosahedral symmetry of genuine icosahedral quasicrystals. However, 
the X-ray images taken in the high-resolution mode do not allow us to locate the examined 
sets. Therefore, as already pointed out in [4], in order to prove practically relevant and rigorous 
results, one has to deal with the non-anchored case of the whole local indistinguishability class 
(or Ll-class, for short) Ll(yl) of a regular, generic icosahedral model set A, rather than dealing 
with the anchored case of a single fixed icosahedral model set A; recall Remark 13.61 for the 
equivalence relation given by local indistinguishability and compare also [18| . 

Remark 4.1. In the crystallographic case of a lattice L in IR^, the Ll-class of L consists of all 
translates of L in IR^, i.e., one has LI(L) = {t+L \ t G R^}. In particular, LI(L) simply consists 
of one translation class. The entire Ll-class LI(ylico(t, VF)) of a regular, generic icosahedral 
model set Aicoit, W) can be shown to consist of all generic icosahedral model sets of the form 
Aicoit,s + W) and all patterns obtained as limits of sequences of generic icosahedral model 
sets of the form Aico(t,s + W) in the local topology (LT). Here, two patterns are e-close if, 
after a translation by a distance of at most e, they agree on a ball of radius 1/e around the 
origin; see [31 [35]. Each such limit is then a subset of some Aico{t,s + W), but s might not 
be in a generic position. Note that the Ll-class Ll(yl) of an icosahedral model set A contains 
uncountably many (more precisely, 2*^°) translation classes; cf. [3] and references therein. 

In view of the complication described above, we must make sure that we deal with finite 
subsets of generic icosahedral model sets of the form ylico(^) s+W), i.e., subsets whose .*-image 
lies in the interior of the window. This restriction to the generic case is the proper analogue of 
the restriction to perfect lattices and their translates in the crystallographic case. Analogous 
to the lattice case [151 [16] and the case of cyclotomic model sets [1], the main algorithmic 
problems of the discrete tomography of icosahedral model sets look as follows. 
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Definition 4.2 (Consistency, Reconstruction, and Uniqueness Problem). Let L = Im(I) 
(resp., L = lo), let W C I^? be a window and let ui, . . . , Um G be m > 2 pairwise non- 
parallel L-directions. The corresponding consistency, reconstruction and uniqueness problems 
are defined as follows. 

Consistency. 

Given functions '■ C^. — > Nq, j G {!,..., m}, whose supports are finite 
and satisfy s\vpY>{Puj) C , decide whether there is a finite set F which is 
contained in an element of T^{W) (resp., T^{W)) and satisfies X^^-F = pu^, 
j G {!,..., m}. 

Reconstruction. 

Given functions Pu^ ■ — > No, j G {1, . . . , m}, whose supports are finite and 
satisfy supp(puj) C , decide whether there exists a finite subset F of an 
element of 2^(VF) (resp., X^{W)) that satisfies Xu^F = puj, j G {1, . . . ,m}, 
and, if so, construct one such F. 

Uniqueness. 

Given a finite subset F of an element of I^{W) (resp., Xj(VF)), decide whether 
there is a different finite set F' that is also a subset of an element of I^{W) 
(resp., Ig{W)) and satisfies XujF = Xu^F', j e {!,... ,m}. 

One has the following tractability result, which was proved for the case of B-type icosahedral 
model sets by combining the results from Section 13.21 with those presented in [1]; cf. [Ml 
Theorem 3.33] for the details. The proof for the F-type case is similar and we prefer to omit 
the straightforward details here. Below, for L G {Im(I),Io}, the L-directions in R H^'^''^'^^ 
will be called L^'^'^'^^ -directions. By Lemma fS. 12( a). the set of Im(I)('^'°'^)-directions and the 
set of Iq^''^'^'* -directions coincide. 

Theorem 4.3. Let L G {Im(I),Io}. When restricted to two L^"^'^'^^ -directions and polyhe- 
dral windows, the problems Consistency, Reconstruction and Uniqueness as defined 
in Definition \4-^ can be solved in polynomial time in the real RAM-model of computation. 

Remark 4.4. For a detailed analysis of the complexities of the above algorithmic problems 
in the B-type case, we refer the reader to [24l Chapter 3]. Note that even in the anchored 
planar lattice case the corresponding problems Consistency, Reconstruction and 
Uniqueness are NF-hard for three or more Z^-directions; cf. [TSl [TB] . 

5. Uniqueness 

5.1. Simple results on determination of finite subsets of icosahedral model sets. 

In this section, we present some uniqueness results which only deal with the anchored case 
of determining finite subsets of a fixed icosahedral model set A by X-rays in arbitrary A- 
directions; cf. Proposition 13. 201 As already explained in Section [U X-rays in non-yl-directions 
are meaningless in practice. Without the restriction to yl-directions, the finite subsets of a 
fixed icosahedral model set A can be determined by one X-ray. In fact, any X-ray in a non- 
yl-direction is suitable for this purpose, since any line in 3-space in a non-yl-direction passes 
through at most one point of A. The next result represents a fundamental source of difficulties 
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in discrete tomography. There exist several versions; compare [2H Theorem 4.3.1], [131 Lemma 
2.3.2], [51 Proposition 4.3], [24l Proposition 2.3 and Remark 2.4] and [23l Proposition 8]. 

Proposition 5.1. Let A be an icosahedral model set with underlying Tj-module L, say A = 
Aico{t-,W^) ■ Further, let U C S'^ be an arbitrary, but fixed finite set of pairwise non-parallel 
L-directions. Then, J-iA) is not determined by the X-rays in the directions of U . 

Proof. We argue by induction on card(C/). The case card(C/) = means U = and is obvious. 
Suppose the assertion to be true whenever card(C/) = k £ No and let card(C/) = k -\- 1. By 
induction hypothesis, there are different elements F and F' of T{A) with the same X-rays in 
the directions of U' , where U' C U satisfies card([/') = k. Let u be the remaining direction 
of U. Choose a non-zero element a G L parallel to u such that a + {F U F') and F L) F' are 
disjoint. Then, F" := {F U {a + F')) - t and F'" := {F' U [a + F)) - t are different elements 
of T[L) with the same X-rays in the directions of U . By Lemma [3.191 there is a homothety 
h: R3 ^ R3 g^pjj ^Yi^^ ^pn ^ pm^ ^ ^^p,f^ ^ ^pm^ ^ ^_ follows that h{F") and h{F"') 

are different elements of .7-"(yl) with the same X-rays in the directions of U; cf. Lemma [2.6[ □ 

Remark 5.2. An analysis of the proof of Proposition ISTTl shows that, for any finite set [/ C 
of k pairwise non-parallel L-directions, there are disjoint elements F and F' of T{A) with 
card(F) = card(F') = 2^=-^) and with the same X-rays in the directions of U. Consider any 
convex subset C of which contains F and F' from above. Then, the subsets Fi := (Cnyl)\F 
and F2 := {C Ci A) \ F' of J^{A) also have the same X-rays in the directions of U . Whereas 
the points in F and F' are widely dispersed over a region, those in Fi and F2 are contiguous 
in a way similar to atoms in a quasicrystal; compare [lH Remark 4.3.2] and [23l Remark 2.4 
and Figure 2.1] (see also [231 Remark 32 and Figure 5]). 

Originally, the proof of the following result is due to Renyi; cf. [32] and compare [2T1 
Theorem 4.3.3]. 

Proposition 5.3. Let A be an icosahedral model set with underlying 'E-module L. Further, let 
U C S'^ be any set of k -\- 1 pairwise non-parallel L-directions, where k £ Nq. Then, J-<,k{A) 
is determined by the X-rays in the directions of U. Moreover, for all F G J^<kiA), one has 
Gfj = F. 

Proof. Let F, F' G !F<k{A) have the same X-rays in the directions of U. Then, one has 
card(F) = card(F') by Lemma (Z^a) and F,F' C by Lemma [M But we have = F 
since the existence of a point in G^\F implies the existence of at least card(C/) > k-\-l points 
in F, a contradiction. It follows that F = F' . □ 

Remark 5.4. In particular, the additional statement of Proposition [5?3l demonstrates that, 
for a fixed icosahedral model set A with underlying Z-module L, the unique reconstruction of 
sets F G T<k{A) from their X-rays in arbitrary sets of k+1 pairwise non-parallel L-directions 
^7 C merely amounts to compute the grids G^. Let A be an icosahedral model set with 
underlying Z-module L. Remark [5.2l and Proposition [5?3] show that T<kiA) can be determined 
by the X-rays in any set of A; -|- 1 pairwise non-parallel L-directions but not by 1 -|- [log2 k\ 
pairwise non-parallel X-rays in L-directions. However, in practice, one is interested in the 
determination of finite sets by X-rays in a small number of directions since after about 3 to 
5 images taken by HRTEM, the object may be damaged or even destroyed by the radiation 
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energy. Observing that the typical atomic structures to be determined comprise about 10^ to 
10^ atoms, one realizes that the last result is not practical at all. 

The following result was proved in [2H Theorem 2.8(a)]; see also [23l Theorem 13(a)]. 

Proposition 5.5. Let d > 2, let R > 0, and let A C be a Delone set of finite local 
complexity. Then, the set 'D^ji{A) is determined by two X-rays in A-directions. 

Since icosahedral model sets A C R'^ are Delone sets of finite local complexity (c/. Re- 
mark [32]), the following corollary follows immediately from Proposition 15.51 in conjunction 
with Proposition 13. 20[ 

Corollary 5.6. Let A be an icosahedral model set with underlying Ti-module L and let R > 0. 
Then, the set T>^fi{A) is determined by two X-rays in L-directions. 

Remark 5.7. Although looking promising at first sight, Corollary 15.61 is of limited use in 
practice because, in general, one cannot guarantee that all the directions which are used yield 
densely occupied lines in icosahedral model sets. 

5.2. Determination of convex subsets of icosahedral model sets. 

Remark 5.8. Proposition 13.201 shows that, for all icosahedral model sets A with underlying 
Z-module -L, the set of L('^''^'^^-directions is precisely the set of yl-directions in Pi H^'^'^'^\ 
Further, by Lemmas 13. 141 and 13. 151 the set of L'^'^-'^'^^-directions maps under ^ bijectively onto 
the set of ^[CsJ-directions. 

The following property is evident. 

Lemma 5.9. Let L G {Im(I),Io}, let U C S'^ be a finite set of L^'^'^'^^ -directions, and let 
F,F' £ J^{t + //('^''^'^)), where t G R^. If F and F' have the same X-rays in the directions of 
U , then $(F — t) and <^{F' — t) have the same X-rays in the directions of ^{U) C S^. 

The following fundamental result follows immediately from [24l Theorem 2.54]; see also [231 
Theorem 15]. 

Theorem 5.10. The following assertions hold: 

(a) There is a set U CS^ of four pairwise non-parallel ZK^]- directions such that, for all 
cyclotomic model sets ylcyc with underlying Tj-module ^[(^5], the set C{Acyc) is deter- 
mined by the X-rays in the directions of U . 

(b) For all cyclotomic model sets /Icyc with underlying 'E-module ^[Cs] and all sets U d'B^ 
of three or less pairwise non-parallel 'L[C,^]- directions, the set C(ylcyc) is not determined 
by the X-rays in the directions of U . 

We are now able to prove the main result of this text by applying the results of pH EH] on 
the determination of convex subsets of cyclotomic model sets with underlying Z-module ^[Cs] 
to the various images ^{{A n (A + H^'^'^'^^)) — A), where A is an icosahedral model set and 
A G yl. 

Remark 5.11. Note that, for a convex subset C of an icosahedral model set A and an element 
A G yl, the intersection C n (A + i?(^'0.i)) is a convex subset of the slice yl n (A + ij{^'0'i)) of 
yl. Hence, $((C n (A + H^^'^'^^)) - A) is a convex subset of $((yl n (A + ii'(^'O'i))) - A). 
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The following fundamental result deals with the anchored case. 

Theorem 5.12. Let L G {Im(I),Io}. The following assertions hold: 

(a) There is a set U C S'^ of four L^'^'^'^^ -directions such that, for all generic icosahedral 
model sets A with underlying Z-module L, the set C{A) is determined by the X-rays 
in the directions of U . 

(b) For all generic icosahedral model sets A with underlying Z-module L and all sets 
{/ C o/ three or less pairwise non-parallel L^'^''^'^^ -directions, the set C{A) is not 
determined by the X-rays in the directions ofU. 

Proof. For part (a), let U' C be a set of four pairwise non-parallel Z[C5]-directions with 
the property that, for all cyclotomic model sets ylcyc with underlying Z-module ^[Cs], the 
set C{Acyc) is determined by the X-rays in the directions of U'. Such a set U' exists by 
Theorem IS.lOr a). We claim that, for all generic icosahedral model sets A with underlying Z- 
module L, the set C{A) is determined by the X-rays in the directions of U := ^'^{U') C 
Cf. Remark 15.81 for the fact that U consists only of L^'^-'^'^^-directions. Assume the existence 
of two different elements, say C and C , of C(yl) having the same X-rays in the directions of 
U. Hence, there is an element Xg A such that Cr\{X + i?(^'0>i)) and C" n (A i^^^'O'^)) are 
different convex subsets of the shce An{X + if(^'O'i)) with the same X-rays in the directions 
of U. By Lemma Ej] and Remark IbTTI it follows that $((C n {X + iJ^^'O'^))) - A) and 
$((Cn(A-Fi?(^'°'^)))-A) are different convex subsets of ^{{An{X + H^'^'^^^^)) - X) having the 
same X-rays in the Z[^5]-directions of U' . Since the set ^>((yln(A+F(^'0'i)))-A) is a cyclotomic 
model set with underlying Z-module l^lCb] by Proposition 13.161 this is a contradiction. 

For assertion (b), let U C be a set of three or less pairwise non-parallel L^'^-'^'^^-directions 
and let yl be a generic icosahedral model set with underlying Z-module L. Consider a slice 
yln(A-Fi?(^'°'^)) of yl, A e A, together with the cyclotomic model set ^>((yln (A-hi?^^'"'^^)) - A) 
with underlying Z-module ^[Cs]; cf. Proposition 13.161 By Theorem 15.101 (b). there are two 
different convex subsets, say C and C , of ^{{A n (A + H^'^'^'^^)) — A) with the same X-rays 
in the Z[C5]-directions of U' := <!>{U) C cf Remark It follows that <!>-^{C) + X and 
$~^(C") -|- A are different convex subsets of (the slice yl n (A -h ij(^'O'i)) of) A with the same 
X-rays in the L(^'°'i)-directions of U. □ 

Remark 5.13. The proof of Theorem 15. 121 shows that the result extends to the set of subsets 
C of generic icosahedral model sets A that are only -convex, the latter meaning that, 

for all A G yl, the sets Cn{X-\- H^'^'^'^^) are convex subsets of the slices yl n (A if(^'°'^)). 

Example 5.14. It was shown in [24l Theorem 2.56 and Example 2.57] (see also [23l Theorem 
16 and Example 3]) that the convex subsets of cyclotomic model sets with underlying Z-module 
^[(■5] are determined by the X-rays in the Z[(^5]-directions of C/5 := {o/|o| | o G 0} C S^, 
where O := {(1 -f r) + Cs, (r - 1) + Cs, -t + C5,2r - Cs} C Z^s] \ {0}. Consequently, as 
was shown in the proof of Theorem 15.121 (a). the convex subsets of generic icosahedral model 
sets yl with underlying Z-module L are determined by the X-rays in the L'-'^'^'^^-directions of 

Remark 5.15. Since, by the work of Pleasants [3l], the Z[C5]-directions of U5 are well suited 
in order to yield dense lines in cyclotomic model sets with underlying Z-module ^[^5], it follows 
that the set of L^'^'^'^^ -directions C/ico from Example 15. 141 is well suited in order to yield dense 
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lines in the corresponding slices AD {X + H^'^'^'^^), A G Tl, of generic icosahedral model sets 
A with underlying Z-module L. In fact, these directions even yield dense lines in icosahedral 
model sets yl as a whole; cf. [31]. In particular, neighbouring slices of the form Ar]{X+H^'^''^'^^), 
A € yl, are densely occupied and hence well separated. Consequently, neighbouring lines in 
any of the directions of U that meet at least one point of a fixed icosahedral model set A are 
sufficiently separated. It follows that, in the practice of quantitative HRTEM, the resolution 
coming from the above directions is likely to be rather high, which makes Theorem 15.121 look 
promising. 

Finally, we want to demonstrate that, in an approximative sense, part (a) of Theorem 15.121 
even holds in the non-anchored case for regular generic icosahedral model sets. Before, we 
need a consequence of Weyl's theory of uniform distribution; cf. [40]. This analytical property 
of regular icosahedral model sets was analyzed in general in [34], [35] and [28]. We need the 
following variant which relates the centroids of images of certain finite subsets of a regular 
icosahedral model set A under the star map to the centroid of its window. 

Theorem 5.16. Let A be a regular icosahedral model set of the form A = ylico(0, VF). Then, 
for all a G R'^, one has the identity 

where A denotes the Lebesgue measure on R'^. 

Proof. This is a consequence of the uniform distribution of the points of A* in the window, 
which gives the integral by Weyl's lemma. The proof of the uniform distribution property for 
model sets can be found in [H [27l [28] . □ 

The following properties of sets U C consisting of four pairwise non-parallel ^[Cs]- 
directions will be of crucial importance: 

(C) For all cyclotomic model sets ylcyc with underlying Z-module ^[Cs], the set C(ylcyc) is 

determined by the X-rays in the directions of U. 
(E) U contains two directions of the form o/\o\, o'/|o'|, where o, a' G ^[Cs] \ {0} satisfy the 

relation 

aoPo'-Poao' G Z[r]^ = {r*|s G Z}, 
where the elements Oq, «o'; /^oi /3o, G are determined hy o = Oq + f^oCb and o' = 
"o' + Po'Cb', cf- Remark [331 

Example 5.17. The set U5 C of four pairwise non-parallel Z[(^5]-directions as defined in 
Example 15 . 141 has property (C) by [24l Example 2.57] (see also [23l Example 3]). Additionally, 
one can easily see that C/5 also has property (E). 

The significance of property (E) is expressed by the following result. 

Proposition 5.18. Let U C 6e a set of four pairwise non-parallel '^[(^5]- directions with 
property (E). Then, for all finite subsets F ofE^^], one has the inclusion 

Proof. This follows from [21 Theorem 1.130] (see also [23l Theorem 12]). □ 
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We are now able to show that, in an approximative sense to be clarified below, for any fixed 
window C whose boundary bd(VF) has Lebesgue measure in R^, the set U^gjB(^4/)C(yl) 

(resp., UyigjP(t^)C(yl)) is determined by the X-rays in any set of L^'^'^^'^^-directions, where 
L = Im(I) (resp., L = Ig), of the form U := ([/'), where U' is a set of four pairwise non- 
parallel Z[C5]-directions with the properties (C) and (E). Since the arguments for the F-type 
case and the B-type case are similar, we present the details for the B-type case only. Let 

F,F' G U C{A), 

say F £ C{Af^^{t,s + W)) and F' G C{Af^^{t' , s' + W)), where t,t',s,s' £ R^, and suppose that 
F and F' have the same X-rays in the directions of U. If F = 0, then, by Lemma [2^ a). 
one also gets F' = 0. One may thus assume, without loss of generality, that F and F' are 
non-empty. Hence, there is an element A G F such that Fn (A-Fif(^'°'^)) and F'n (A-h-ff^^'^'^^) 
are non-empty finite sets with the same X-rays in the directions of U. Then, by Lemma \5M 
the non-empty finite subset $((Fn(A-hij(^'°'^)))-A) ofZ[C5] (c/. LemmaEH]) and the non- 
empty finite subset $((F' n{\ + - A) of C have the same X-rays in the four pairwise 
non-parallel Z[(^5]-directions of ^{U) = U' . Then, by Lemma [T4l and Proposition 15.181 in 
conjunction with property (E), one obtains 

<I>((F n (A + F(-'O'i))) - A), <I>((F' n (A + - A) C g'J«^"(^+^*""''))-") c ZKs] • 

Thus, one gets 

(6) Fn(A + //("'0'^)),F'n(A + F("'0'^)) c t + L. 

Since F' n (A + ct' + L, Relation {6]) implies that t + L meets t' + L, the latter being 

equivalent to the identity t + L = t' + L. Note also that the identity t-|-L = i'-|-L is equivalent 
to the relation t' — t £ L. Clearly, one has 

F-t £ C{Af,,iO,s + W)). 

Moreover, since the equality 

Atit' -t,s' + W)= Al,iO, is' + it' - ty) + W) 

holds, one further obtains 

F'-t £ CiAl,{t'-t,s' + W)) = CiAl,{0,{s' + it'-tr) + W)). 

Clearly, F — t and F' — t again have the same X-rays in the directions of U . Hence, by 
Lemma [2?2T b). F — t and F' — t have the same centroid. Since the star map .* is (Q-linear, it 
follows that the finite subsets {F — t)* and {F' — t)* of R'^ also have the same centroid. Now, 
if one has 

F-t = Bnia) n ylL(0,s + Ty) 

and 

F' -t = BR,{a') n Af,,iO,{s' + it' -ty) + W) 

for suitable a, a' £ R'^ and large R,R'>0 (which is rather natural in practice), then Theo- 
rem 15.161 allows us to write 
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voWy,+,^ ^^^^^^^ ^ card (F - t) 



X 



\ V X* 

card iF' - t) ^ 

1 



vol(VF) 



2/dA(y) . 

{s>+{t'-tY)+w 



Consequently, 



/ ydA(y) « {s' + {t'-tr)+ [ ydX{y). 
Jw Jw 



and hence s s' + [t' — t)* . The latter means that, approximately, both F — t and F' — t 
are elements of the set C{A^^^{Q,s + W)). Now, it follows in this approximative sense from 
property (C) and Theorem EE] that F -t^ F' - t, and, finally, F w F' . 

Remark 5.19. The above analysis suggests that, for all fixed windows W C I^? whose 
boundary bd(W) has Lebesgue measure in R^, the sets of the form ^A&xf{w)^{-^) (resp., 
U^gjF(j^)C(yl)) are approximately determined by the X-rays in the four prescribed L^"^'*^'^)- 
directions of C/ico, where L = Im(I) (resp., L = Iq); c/. Examples 6.12 and 6.15. Additionally, 
in the practice of quantitative HRTEM, the resolution coming from the directions of U\co is 
likely to be rather high, which makes this approximative result look even more promising in 
view of real applications; cf. Remark 15.151 



6. Outlook 

For a more extensive account of both uniqueness and computational complexity results in 
the discrete tomography of Delone sets with long-range order, we refer the reader to [24]. This 
reference also contains results on the interactive concept of successive determination of finite 
sets by X-rays and further extensions of settings and results that are beyond our scope here; 
compare also [23]. Although the results of this text and of [24] give satisfying answers to the 
basic problems of discrete tomography of icosahedral model sets, there is still a lot to do to 
create a tool that is as satisfactory for the application in materials science as is computerized 
tomography in its medical or other applications. First, we believe that it is an interesting 
problem to characterize the sets of yl-directions in general position having the property that, 
for all icosahedral model sets A, the set of convex subsets of A is determined by the X-rays in 
these directions; compare [131 Problems 2.1 and 2.3]. Secondly, it would be interesting to have 
experimental tests in order to see how well the above results work in practice. Since there is 
always some noise involved when physical measurements are taken, the latter also requires the 
ability to work with imprecise data. For this, it is necessary to study stability and instability 
results in the discrete tomography of icosahedral model sets in the future; cf. [1]. 
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